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BODY ROTATION EFFECTS ON MELTING ABLATION 


by Thomas H. Cochran and Simon Ostrach* 

Lewis Research Center 

SUMMARY 

The effects of rotation on the ablation of an axisymmetric body composed of a highly 
viscous material and subject to deceleration in an airstream were studied analytically 
and experimentally. The tests were conducted so that a parameter representing the im- 
portance of rotation effects and which was obtained from a nondimensionalization of the 
governing equations was varied from 0 to 4. 2. 

The results indicate that, although rotation had no effect on the ablation process in 
the vicinity of the stagnation point of the body, significant changes in the character of the 
generated waves and in the flange were caused by rotation. Liquid was thrown from the 
body in the vicinity of the flange when the spin parameter was 0. 88 or greater. 


INTRODUCTION 

With the advent of space travel, the need has arisen for the protection of payloads 
from the high heat fluxes encountered during reentry into the Earth T s Atmosphere. Solutions 
to this problem have centered around the use of external shields which are sacrificed in 
an ablation process in order to protect the vehicle interior. For ballistic applications 
which have relatively short-time high -heat -intensity reentry characteristics, ablative 
materials that sublime have been of greatest interest. However, for the longer time, 
lower heat intensity reentry environments which a manned spacecraft may encounter, 
ablative materials can melt before they vaporize. Therefore, numerous research pro- 
grams, which are discussed in the following section, have been initiated to study the fluid 
dynamics and heat transfer of thin melt layers. 


* Professor of Engineering, Head of Division of Fluid, Thermal, and Aerospace 
Sciences, Case Institute of Technology, Cleveland, Ohio. 




Figure l. - External structure of tektites. 
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Another area in which melting ablation has become of significance is in the explan- 
ation of tektites. These small (about 2 cm in diam) button-shaped glassy bodies (see 
fig. 1) have been found over large areas at numerous sites around the world and bear no 
geological resemblance to their surroundings. On the surface of some of these objects, 
such as those found in Australia, are ring-shaped waves and an equatorial flange. 
Chapman (ref. 1), and Ostrach (ref. 2) have found that this distinctive appearance can 
be produced on bodies ablated in a wind tunnel. Based on such aerodynamic tests and 
exhaustive geometric and structural investigations, Chapman has proposed that tektites 
are extraterrestrial, having ablated during entry into the Earth's atmosphere. Further, 
by measuring the amount of ablation at the stagnation point and the spacing of the ring 
waves, he has obtained trajectory information that has led him to confirm the suggestion 
that the origin of tektites is the moon. O'Keefe (refs. 3 and 4) agrees with a lunar 
origin. However, he believes that the individual tektites did not come directly from the 
moon as Chapman has proposed but were torn from a larger body as it passed through 
the Earth's atmosphere. Regardless of their nature before reaching the Earth, tektites 
provide a convenient model and standard for melting ablation research. 

A possible condition imposed on a body entering the atmosphere is that of spinning. 
Chapman (ref. 1) suggests that tektites rotate during their ablation; however, no mention 
is made of expected rates or effects on the ablation process. It is possible that rotation 
effects (centrifugal and coriolis forces) could play an important role in the dynamics of 
a melt layer, and hence the ablation process. This suggests that the variation in the 
amount of ablation apparent on different tektites may be attributed to rotation rather than 
trajectory characteristics. 

This report presents the results of an investigation of the effects on melting ablation 
of rotation about the major axis of a body of revolution similar to that of a tektite. The 
problem was, first, treated analytically to obtain the pertinent governing equations and 
then experimentally to measure the effects of rotation on the ablation process. Test 
models were spun over a wide range of speeds (0 to 400 rpm) in an effort to obtain a 
clear picture of rotation effects on the melt. 


PREVIOUS ABLATION WORK 

The early work in ablation was keyed to design criteria because of the need for 
immediate information. Approximate analysis and experimentation were used to provide 
the engineer with the necessary tools to treat the problem grossly. Adams has provided 
a survey of these design techniques in reference 5. 
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Early attempts to treat aerodynamic ablation in a more exact manner were confined 
to the stagnation region. Sutton's work in reference 6 is considered a classic of these 
initial efforts because he was able to show the importance of the melt layer in the ablation 
process. However, he compromised his solution by neglecting transient effects. 

A significant advance in defining melting ablation phenomena was achieved by 
Ostrach, Goldstein, and Hamman in reference 7. Here, the problem was approached in 
a general form by including transient effects and vehicle deceleration. Their solution, 
which was carried out around the body as well as at the stagnation point, showed the 
importance of body forces on the melt and predicted a wave-like behavior. However, 
only a single wave was indicated, rather than the repeated waves found on the tektites. 
Also, the solution was essentially quasi-steady. 

McConnell (ref. 8) extended Ostrach’s work by applying a perturbation analysis to 
the early stages of ablation where conduction dominates convection. He also improved 
Ostrach's quasi-steady solution by generalizing the boundary condition. Comparison of 
his theoretical results with experimental data taken with asphalt spheres in a wind tunnel 
indicated good agreement. 

In recent work in this problem, Chen (ref. 9) has extended the detailed analysis of 
the liquid-layer equations by Ostrach and McConnell to include the gas -boundary -layer 
equations. By considering the coupling of the gas and liquid equations at their interface, 
he has been able to obtain pertinent reference quantities which permit a more rigorous 
nondimensional treatment of the governing equations. This approach to the problem also 
improves the validity of boundary conditions at the interface. 

The value of Chen's analysis has been demonstrated to some extent by the work of 
Edling (ref. 10). He investigated the characteristics of waves generated on an ablating 
flat plate for several different materials. A correlation of the results was accomplished 
through an ablation parameter which was derived by Chen and predicted earlier by 
Ostrach and McConnell (ref. 2). This parameter relates the rate at which a substance 
melts to the incident heating rate. 


ANALYSIS 

Melting Ablation Model 

The mathematical model assumed is that of an axisymmetric body composed of a 
material which, when in a liquid state, has a viscosity that is highly temperature- 
dependent. The thermal conductivity, specific heat, and density of the material will be 
considered to be independent of temperature. The body is assumed to decelerate at a 
constant rate in the atmosphere and is subject to heating sufficient to cause a layer of 
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melted material to flow along the surface. Because of the viscous nature of the material, 
the layer is very thin compared with the body size, and no distinct liquid-solid interface 
is identifiable. The flow is driven by aerodynamic, gravity, and rotational forces. 
Evaporation of material is neglected. 

A schematic drawing of the model is shown in figure 2. Since the liquid layer is thin 
compared with the radius of curvature of the body, a Cartesian coordinate system is 
chosen normal to and along the liquid-gas interface. The x-axis is along the interface 
meridionally, the y-axis is perpendicular to the interface, and the z-axis is along the 
interface azimuthally. The coordinate system is fixed at the liquid-gas interface and is 
rotating at a constant rate equal to that of the solid portion of the body. It is assumed 



(a) Three-dimensional view. 
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that acceleration terms associated with the unsteady motion of the interface relative to 
the body are negligible. 


Governing Liquid Equations 

Dimensional equations. - The equations governing the motion of the liquid layer 
Continuity: 


JL (RU) + JL (RV) = 0 
ax ay 


Momentum: 


9U + U 9U + V^-A 


at 


ax 


3Y 



ap a 

ax 3Y 




'av + auV 
v ax dY ) 



Pi 







ap + _a_ 
aY ax 







Pi 


9W + U ^ + V^ + 2Uco^ 

at ax aY dx 


2Vu> 



Energy (neglecting thermal expansion): 


M C p> 


/9T + U ar + V ar' 
l at ax SY ' 




+ $ 


are: 

( 1 ) 


( 2 ) 


(3) 


(4) 


(5) 
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where <i> is the dissipation function and is expressed as 


* + H 


au 

9X 


+ 2 


dv 

9Y 


+ 2^u av + 

9Y 9X 


9U\ 2 + /9V\ 2 + /9W\ 2 + /9W x2 
9Y/ \9X/ \9Y / \9X 


( 6 ) 


Constitutive relation: 


T - T 


,-N 


Hi 


( 7 ) 


(All symbols are defined in appendix A. ) 

In writing these equations, the tangential velocity due to spinning is taken to be 
independent of y. The error involved in assuming this is negligible because of the 
thinness of the liquid layer. Deceleration, coriolis, and centrifugal terms arise in the 
transformation from a stationary coordinate system to the one fixed at the liquid-gas 
interface. 

Nondimensionalization . - Simplification of the equations is accomplished by a nondi- 
mensional technique. By careful selection of characteristic values for the variables, the 
relative importance of the different terms in the equations can be assessed. The nondi- 
mensional variables are as follows: 


X Y U V 

x = — , y = — , m = , v = 

L 6jL U 


ref A l U ref 


w = 


W_ 

A 2 U ref 


H 


R 


O I t JL Vr 

n = f t =-, r = — , y T = 


i l 


T 7T’ 
e Z L 


p = 


P U 2 

^cr _oo 


T - T 


e 


T i - 

1 °0 


Of the characteristic values, the scaling parameter for the liquid velocity layer thick- 
ness 6^, the reference velocity U rg j, the characteristic time value a, the scaling par- 
ameter for the liquid thermal layer thickness e^ , the scaling parameter for the velocity 
normal to the liquid-gas interface A^, the temperature at the liquid -gas interface T^, 
and the scaling parameter for the velocity in the azimuthal direction Ag are unknown 
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and must be determined from an order -of -magnitude analysis of the nondimensional 
equations and boundary conditions . Chen (ref. 9) evaluated 6^ , , T., and Aj, 

and his analysis is presented in detail in appendix B. The remaining parameter Ag 
must be obtained from the equation of motion for the azimuthal direction. This equation 
in nondimensional form is 



Assuming that any motion of the liquid in the aximuthal direction relative to the spinning 
reference system is driven by coriolis forces and opposed by viscous forces results in 



After inserting the new variables and gathering terms, the equations are 

— (ru) + — (rv) = 0 (10) 

3x dy 
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A list of the equation parameters and their values for a possible test material at a 
condition of interest is presented in table I. Physical properties used in computing values 
for table I are given in table II. Examination of the parameters reveals some which are 
not in common usage in the scientific community, and, therefore, warrant description: 

N ratio of pressure driving effects to viscous retarding effects, will be referred to 

r 

as pressure parameter 

N g may be interpreted as the ratio of centrifugal inertia effects to viscous effects, and , 
essentially, is a modified form of a Taylor number 

g modified form of reciprocal Froude number 

E modified form of Eckert number 

The ratio N g /Np, therefore, relates the importance of centrifugal driving effects to 
pressure driving effects in the liquid layer. 
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TABLE I. - TYPICAL PARAMETERS FOR 650 OIL AND TEKTITES 


Parameters 

Material 


650 Oil a 

Tektite 

Liquid Reynolds number, Re^ 

3. 04xl0' 3 

3.82X10' 3 

Rotation Reynolds number, Re^ 

1. 97x10“ 1 a) 

0. 42X10' 2 w 

Liquid scaling parameter for 

0. 5x10“ 1 

0.71X10' 1 

direction perpendicular to 



liquid-vapor interface, 6^ 



Modified Taylor number, N g 

3. 19X10' 2 w 2 

2.3xl0‘ 5 w 2 

Dimensionless pressure 

7.5 

1.18 

parameter, 



Liquid thermal layer thickness, 
e. 

0.665X10' 1 

1.27X10' 1 

Liquid Prandtl number, Pr^ 

3.95X10 4 

2.85X10 4 

Modified Eckert number, E 

1. 08xl0" 6 

7.35X10' 6 

Dimensionless acceleration 

-1. 19 

-0. 43 

parameter, g 




a Values are based on a model of 10" 1 ft (3. 05 cm) radius in a 
wind tunnel for an air velocity of 45 ft/ sec (1. 37 m/ sec). 
The tektite calculations are based on estimated entry 
conditions as outlined in appendix B of ref. 8. 


Order -of -magni tu de analysis . - An order -of -magnitude analysis based on the values 
for the parameters contained in table I results in the following simplified equations: 


— (ru) + — (rv) = 0 (15) 

9x 3y 
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TABLE n. - PROPERTIES OF 650 OIL, AIR, AND TEKTITES 


Property 


Substance 



650 Oil 

Air 

Tektites 

Temperature properties evaluated at- 





°F 

60 

78 

5300 

3040 

°K 

289 

299 

3203 

1943 

Thermal conductivity, k: 





Btu/(ft)(sec)(°F) 

3.3X10 -5 

4. 16X10" 6 

9.25X10" 5 

2. 57X10" 4 

j/(m)(sec)(°K) 

2. 05x10" 1 

2. 59xl0“ 2 

5.75X10' 1 

1.6 

Specific heat, Cp: 





Btu/(slug)(°F) 

14.9 

7.7 

36.8 

7.35 

j/(kg)(°K) 

5. 95xl0 2 

3. lxlO 2 

1. 48xl0 3 

2.96X10 2 

Density, p: 



10" 6 


slug/ft 3 

1.73 

2.3X10' 3 

4.65 

kg/m 3 

8.9X10 2 

1. 185 

5. 15X10" 4 

2.4X10 3 

Viscosity, Jl: 


0.39X10" 6 

1. 61xl0" 6 


slug/ (ft) (sec) 

0.88X10" 1 

1.00 

kg/(m)(sec) 

4.22 

1.87X10' 5 

7.72X10' 5 

4. 79X10 1 

Melting temperature: 
°F 

15 

264 



1500 

1089 

°K 



Source 

Ref. 9 

Standard 

conditions 

Ref. 11 

Ref. 8 


9p 

3y 




(17) 


u 


dr 

dx 



(18) 


de de de / i\ d 2 e 

at ax ay \p) dy 2 


(19) 
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( 20 ) 


M 


l 


= e 


-N 


o 

where /3 = Pr^ Re^ and relates the importance of convection effects to conduction 
effects. 

These equations are the same as those obtained by Chen with the addition of centrif- 
ugal inertia terms in the equations of motion. The presence of these terms is illustrated 
in figure 2 by a force diagram. It is evident that in the meridional direction X, the 
centrifugal force aids in driving the liquid along the surface. For a spherical body, the 
effect would be negligible at the stagnation point because of the dependence of the centrif- 
ugal force on distance R from the rotation axis and at the equator because of the geom- 
etry of the body dR/dX = 0. The maximum driving effect would occur at a position 
equidistant from these two minimum points. It should be noted that the centrifugal force, 
which is also a body force, acts in the opposite meridional direction compared with the 
deceleration force. Therefore, for high rotation rates, at some meridional positions, 
the centrifugal force may become sufficiently large to offset or overpower the deceler- 
ation force and its effects on the dynamics of the melt layer. 


Governing Gas Equations 

The problem of an axial airstream flowing over a rotating axisymmetric body has 
been treated by Schlichting in reference 12. An analytical investigation yielded boundary- 
layer parameters, frictional drag, and torques as a function of the ratio of circumfer- 
ential velocity to free-stream velocity. The results for the case of a spherical body 
indicate that, for spin rates pertinent to this work, the momentum loss thickness, torque, 
and frictional drag are negligibly affected by rotation. Effects of rotation on the heat 
transfer from a body of revolution in an axial airstream has been studied experimentally 
by Lewis and Ruggeri (ref. 13). Their results indicate that, for tangential velocities up 
to 105 feet per second (32 m/sec) and free-stream velocities of approximately 250 miles 
per hour (112 m/sec), no significant change in the heat-transfer coefficients occurred. 

It is evident, therefore, that, for the problem at hand, the governing equations for 
the gas may be approximated with little error by those for air flowing over a nonrotating 
body. This formulation has been accomplished by Chen whose equations, in ordered form, 
are presented here for the sake of completeness: 

— (ru) + — (rv) = 0 (21) 

9x 0y 
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( 22 ) 


u J“ + v^ = . 9P + 

3x 9y 9x gy2 

i?E = 0 

sy 

u^ + v ^= / 1 V 2g 

9x 3y \Pr J ~ 2 

J \ g/ 9y 


Tektite Parameters 

Because the ordering procedures were carried out for a specific test material and 
conditions, it is of interest to compare the calculated parameters with those obtained for 
a tektite entering the atmosphere. This has been done in table I for estimated entry con- 
ditions as presented in reference 8. Physical properties of air and tektites for these con- 
ditions are shown in table n. Comparison of the terms in table I indicates that, with the 
exception of the rotation parameters Re^ and N g , the values of a parameter for 650 oil 
and tektites are of the same order of magnitude. 


(23) 


(24) 


EXPERIMENTAL APPARATUS 
Test Model 

The characteristics deemed desirable in the material composing the model were 

(1) That its viscosity be highly dependent on temperature 

(2) That it have a relatively high specific heat and low thermal conductivity 

(3) That its melting point be relatively low 

This last characteristic was particularly advisable in that it removed the necessity for a 
high-temperature aerodynamic facility. The material chosen was 650 oil, the properties 
of which were given earlier in table n. 

The oil was formed on the exterior of a 1. 1-inch- (2. 8-cm-) diameter wooden sphere. 
Attached to the sphere was a 0. 375-inch- (0. 95-cm-) diameter, 3-inch- (7. 62-cm-) long 
plastic tube. These materials were chosen to minimize conduction away from the oil 
through its supporting structure. The method used to form the oil on the wooden sphere 
was developed by Chen (ref. 9). By alternately dipping the sphere in oil and liquid 
nitrogen, thin layers of solidified oil were built up until the desired size was reached 
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(approximately 2 in. (5. 08 cm) in diam). Figure 3 schematically shows the details of the 
model. 

Temperature within the model was measured by a copper -constantan thermocouple 
located adjacent to the wooden sphere and on the major axis of the model. 


Aerodynamic Facility 

To simulate deceleration through the atmosphere, the experiments were conducted 
in a vertical air duct. This wind tunnel was designed by Chen and is shown in figure 4. 
A centrifugal blower, driven by a 7. 5-horsepower (5580-W) induction motor, forced 
room -temperature air through the tunnel at approximately 45 feet per second (13. 7 
m/sec). As illustrated in the schematic drawing, guide vanes, baffles, and screens 
were used to reduce the turbulence level in the flow and provide a flat velocity profile in 
the vertical portion of the tunnel. The ducting was constructed of sheet metal with the 
exception of one section which was clear plastic so that the ablating model could be 
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Figure 4. - Aerodynamic facility. 
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, observed and photographed. Motion pictures were taken with a 16-millimeter camera 
that was electrically driven at 25 frames per second. Illumination was supplied by two 
photoflood lights. Time reference was provided by a digital clock mounted outside the 
tunnel in the field of view of the camera. 

Apparatus to spin the models was installed in the tunnel as shown in figure 5. A 
1/ 15-horsepower (0. 496-W) variable -speed direct-current motor provided the driving 
force through a system of gearing. At the motor end of the power train, worm gearing 
was used for the slower rotation speeds, while a system of spur gears was used to 
attain faster rates. The speed of the motor was checked under test conditions with a 
stroboscope and varied a maximum of 1^- percent. 
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The test models were fixed on the spinning apparatus by slipping the plastic tube 
over the end of the vertical shaft and shoving upward until a firm friction fit was effected. 
Leads from the thermocouple were connected to two slip rings, one copper and one 
constantan. Brushes rode on the rings and enabled the model temperature to be recorded 
on a paper chart recorder which operated at a constant speed. The cold junction was 
maintained at 32° F (273° K). 


Preliminary Testing 

After the apparatus to spin the models was installed in the tunnel, velocity profiles 
were obtained in the A -A and B-B directions, which are shown in figure 6. This was 
done in planes immediately above and below the test model location to ensure that the air- 
flow was not being adversely affected in the vicinity of the model. Measurements were 
made with a pitot tube which could traverse 10 inches (25. 4 cm) into the tunnel. Pres- 
sure was recorded on a micromanometer that was "zeroed" at barometric conditions. 

The results are plotted in figure 6. It should be noted in this figure that the grid assumes 
a double role: (1) distance and (2) percent deviation of the velocity from that measured 
at a particular position. This position was selected in the top plane of measurement on 
what was the major axis of the model when it was installed on the spinning apparatus. 

The curves show that, within a distance of 2 inches (4. 08 cm) from the major axis of the 
model, the air velocity varied a maximum of ±2 percent. 

Turbulence intensity at the center of the test section was measured with a hot wire 
anemometer and found to be ±0. 1 percent. Air temperature in the tunnel was measured 
with a copper -constantan thermocouple and recorded on a chart recorder. 


TEST PROCEDURES AND DATA ACQUISITION 
Procedures 

After a model had been made, it was placed in a constant temperature bath. The 
bath consisted of a Dewar half filled with liquid nitrogen and covered with a heavy cloth. 
The model was suspended in the top half of the Dewar from a laboratory stand, and the 
thermocouple leads were attached so that the temperature in the model could be 
monitored on the chart recorder. After 5 to 10 minutes in the bath, the temperature de- 
crease became slow (approximately 0. 5 F^min (0. 28 K°/min)). In his experimental work, 
Chen (ref. 9) used a model with thermocouples at various positions in the oil and found 
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that fairly uniform conditions were attained after 5 minutes in the bath. When these 
approximate equilibrium conditions were attained, 

(1) The model was removed from the bath and placed in the wind tunnel. 

(2) The spinning rate was set at the desired speed and the blower started. 

(3) Motion pictures were taken intermittently during the test which generally required 
about 15 minutes. Time scales were checked and correlated by occasionally recording 
the digital timer reading on the chart recorder. 

(4) When all the oil had been ablated from the wooden sphere, the model was removed 
from the tunnel, and pressure measurements were made with a pitot tube. 


Data Reduction 

The zero point on the time scale was chosen as that time when the model was fixed 
on the spinning apparatus in the wind tunnel. This was exactly determined from the 
chart recorder because the model thermocouple started to ’’read out” at this point. 

Analysis of the film was accomplished on a motion-picture analyzer that magnified 
the image eight times. Step by step film advance and vernier scales permitted accurate 
measurements to be made. In addition to observing the general ablation process, 
recession of the stagnation point and wave crest motion were determined. Also, profiles 
of the models were obtained by projecting the film onto paper and tracing the basic 
model outline (waves were not included). 

Calculation of the pertinent nondimensional parameters for each test required, first, 
that the characteristic values for the variables be determined. Initially, the interface 
temperature T. and the liquid viscosity at the interface were obtained by a trial - 
and-error solution of equations (B21) and (B22). This involved selecting a value of T^ 
from which a value of could be obtained from the viscosity -temperature data of 

Edling (ref. 10). The value for was then inserted into equations (B24) and (B25), 
and a value for T^ was calculated. The procedure was continued until the assumed and 
calculated values of T^ were the same. The remainder of the characteristic values 
could then be obtained from their formulation as presented in appendix B. 


RESULTS AND DISCUSSION 
Test Conditions and Parameters 


A summary of the test condtions and equation parameters for the experimental runs 
is presented in table IH. It can be seen that all the major parameters in the momentum 
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TABLE m. - TEST CONDITIONS AND PARAMETERS 


[Gas Prandtl number, 0. 72. ] 


| Condition or parameter 



Experimental run 


| 


1 

2 

3 

4 

5 

6 

7 

Rotation rate, rad/ sec 

C 

1.51 

8.2 

13. 

19.4 

3C 

42.6 

Gas free-stream temperature, T : 








°F 

77.8 

77.8 

77.8 

76.8 

76.8 

76.8 

76.8 

°K 

298.6 

298.6 

298.6 

298. C 

298.0 

298.0 

298. 0 

Gas free-stream velocity in X-direction, 
IT 








-<K> * 

ft/ sec 

46 

45.6 

45.6 

45.2 

45.2 

45.2 

44.7 

m/sec 

1. 40 

1.39 

1.39 

1.38 

1.38 

1.38 

1.36 

Center -of -model temperature, T^: 








°F 

-1.5 

-2.6 

-0.2 

-4.0 

-5.2 

-0.2 

-0.2 

°K 

254.5 

253.9 

255.3 

253.2 

253. 1 

255.3 

255. 3 

Length characteristic of model size, L: 








ft 

0. 082 

0.082 

0. 076 

0. 076 

0. 079 

0.078 

0.078 

m 

0. 025 

0.025 

0. 023 

0. 023 

0.024 

0. 024 

0. 024 

Gas Reynolds number, Re^ 

2. 10X10 4 

2 . 21 X 10 4 

2. 04xl0 4 

2. 02X10 4 

2. 10X10 4 

2 . osxio 4 

2. 06x1 0 4 

Liquid Prandtl number, Pr^ 

1. 26x10® 

1. 34xl0 5 

1. 18x10® 

1. 42x10® 

1. 51x10® 

1. 34x10® 

1. 34xl0 5 

Ablation parameter, y 

0.530 

0. 526 

0.542 

0.510 

0. 500 

0. 526 

0. 526 

Liquid interface viscosity, /a ^ ; 








slug/ (ft)(sec) 

0.282 

0.298 

0. 263 

0.316 

0.334 

0.298 

0. 298 

kg/(m)(sec) 

13.5 

14.3 

12.6 

15.2 

16.0 

14.3 

14.3 

Interface temperature, T. : 








°F 1 

50.25 

50. 00 

50. 50 

49.75 

49, 50 

50.00 

50. 00 

°K 

283.3 

283.2 

283.4 

283.0 

282.9 

283.2 

283.2 

Reference velocity in X-direction, U f : 








ft/ sec 

0.76X10" 3 

0.70X10" 3 

0. 82xl0' 3 

0.67X10" 3 

0. 62X10" 3 

0.72X10' 3 

0.71X10" 3 

m/sec 

2.36xl0‘ 4 

2.14X10' 4 : 

2.50X10' 4 

2. 04x10“ 4 

1. 89xl0“ 4 : 

2.20X10' 4 

2. 16X10" 4 

Liquid scaling parameter for direction 

0. 0827 

0. 0794 

0. 0854 

0. 0848 

0. 0817 

0. 0846 

0. 0852 

perpendicular to liquid-vapor interface, 
6 7 








l 

Time characteristic value, <j 

1. 02X10 2 

1.17X10 2 

0. 93X10 2 

1. 13X10 2 

1.28X10 2 

1. 08xl0 2 

1. 10X10 2 

Convection parameter, j3 

0.325 

0.280 

0. 354 

0.284 

0.254 

0.312 

0.315 

Absolute value of dimensionless acceler- 

0.93 

0.94 

0.88 

0.90 

0.93 

0.92 

0.94 

ation parameter, g 








Dimensionless pressure parameter, N p 

12.0 

11.8 

12.2 

12.0 

11.9 

12.2 

12.2 

Dimensionless spinning parameter, N g 

0 

0. 066 

1.74 

4.7 

10. 5 

24.6 

51.0 

Spin parameter, N g /N p 

0 

0. 001 

0. 14 

0. 40 

0. 88 

2.20 

4.20 
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equations were kept relatively constant except the Taylor number N , thus permitting a 

s 

clear picture of rotation effects to be obtained. The spin parameter N /N ranged 

s p 

from 0 to 4. 2 (u> = 0 to 400 rpm) while the pressure parameter Np and reciprocal 
Froude number g were maintained at approximately 12 and -0. 9, respectively. 


General Observations 

Wave motion . - The characteristics and motion of the generated waves were closely 
observed on the film analyzer. In the case of the stationary model, the waves initially 
formed concentrically about the stagnation point. As these waves moved back along the 
model surface, they generally broke up, with some portions of the waves moving faster 
than others. At large displacements from the stagnation point, the waves were, in 
general, highly distorted, having lost their axisymmetric nature entirely. This is shown 
for typical waves in the photographs in figure 7(a). 

For the rotating models, the waves again initially formed concentrically about the 
stagnation point and broke up as they moved back along the surface. However, it became 
evident that, as the rotation rate was increased and at large meridional X displacements, 
the waves were more axisymmetric. As shown in figure 7(b), the waves had the shape of 
a helix of one revolution. It is apparent, therefore, that the net effect of rotation was to 
make the waves more symmetric about the major axis of the model. 

A possible cause for wave distortion was the imperfect velocity profile (±2 percent 
deviation from the stagnation point). To ascertain this effect, a test was conducted at a 
low rotation rate (14 rpm). For such a condition, the centrifugal effects would be small 
(N g /Np = 0. 001), and the small discrepancies in the velocity profile should be evened out 
about the model. Observation of the data film revealed no significant change in the waves 
from those observed on the stationary model. 

Wav e amplitude . - Another change in the wave characteristics which was observed in 
the data films (see fig. 7) concerned the wave amplitude. It was apparent that the 
amplitude of the waves decreased as the. rotation rate was increased. However, attempts 
to measure the reduction were fruitless because of the minuteness of the amplitude itself. 
An explanation for this reduction is concerned with the opposing centrifugal and decelera- 
tion forces in the melt layer, as has been discussed. Ostrach (ref. 14) indicates that the 
presence of a deceleration force in the melt layer results in the amplitude of the waves 
becoming greater. Therefore, by rotating the models and negating the deceleration 
effects with the centrifugal effects, it seems reasonable that the wave amplitude should 
decrease. 
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(a) Spin parameter, 0. 



(b) Spin parameter, 4.2. 


Figure 7. - Appearance of models during ablation. 






Recession of Stagnation Point 


The progress of the ablation process in the vicinity of the stagnation point may be 
determined by the recession of that point from its original position. This measurement 
was made for each of the test runs except run 3. During this test, the digital clock 
stopped for a period of time that could not be accounted for. 

The results are presented in dimensional form in figure 8 and nondimensionally in 
figure 9. The shapes of the curves agree with the experimental and analytical results 



Figure 8. - Recession of stagnation point for various spin rates. 



Dimensionless time, t 

Figure 9. - Dimensionless recession of stagnation point. 
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of reference 8, that is, the ablation velocities, as indicated by the slopes of the curves, 
are small initially and approach a steady-state value after some time has elapsed. It is 
apparent from figure 9 that, with the exception of data from the stationary model, Chen's 
nondimensional technique correlated the results. The discrepancy in the results for the 
one test may well be accounted for by experimental error. Primarily, the conclusion to 
be drawn from these data is that no clear effect of rotation was noted in the recession of 
the stagnation point. This would be expected because of the dependence of centrifugal 
force on distance from the axis of rotation. 
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Body Shape Characteristics 


Front surface . - In order to study further the effect of rotation on the ablation pro- 
cess in the vicinity of the stagnation point, the profiles of the models at various times 
were investigated, as indicated in Data Reduction. The results are presented in figure 10 
for a cross section of the test conditions. What are considered to be the "front" sur- 
faces of the models are denoted in the figure. This area was defined as being upstream 
of the meridional X inception point of the waves as generated on the stationary model. 
The need for this specification is brought out in a later discussion (see section, Wave 
Dynamics). Investigation of these front areas reveals that the surface shapes remained 
essentially spherical during the ablation process for all test conditions. Apparently, for 
the meridional displacements concerned, the centrifugal effects do not become large 
enough to affect the ablation process. 

Flange . - A significant change that occurred in the body shape as a function of 

rotation rate was observed in the vicinity of the flange. It is shown in figures 7 and 10 

that, for no rotation and small rates of rotation, the melted material accumulated in an 

unsymmetrical mass in the vicinity of the equator of the models. Occasionally, relatively 

large chunks of liquid were torn from this flange by the shearing effect of the airstream. 

At higher rotation rates (N /N > 0. 88), the flange was swept back along the model 

s p 

surface because of the centrifugal effects acting normal to and along the liquid-vapor 
interface. In addition, liquid was propelled from the model in a small steady stream 
around the entirety of the flange. 


Wave Dynamics 

The effects of rotation on the ablation process in the vicinity of the stagnation region 

and flange having been investigated, attention was focused on the surface area between 

these positions by studying the motion of the generated waves. Consequently, the crest 

position of the waves as a function of time was measured for the different test conditions. 

The data for the test in which the spin parameter N /N equaled 0. 001 were used 

s p 

instead of that from the stationary test because the entire circumference of the waves 
could be viewed in the former and not in the latter. Since the waves were not exactly 
axisymmetric, a maximum and minimum position at any one time was measured so that 
an average meridional X displacement could be determined. Waves were selected from 
that portion of the film when the models were ablating approximately at a constant rate, 
as determined from the stagnation point recession data (t > 3. 5). Another prerequisite 
in selecting waves was that their entire lifetime from inception until they merged with 
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the flange be visible. Figures 11(a) to (f) present the results for the 28 waves measured. 
In these plots, each curve represents a different wave. 

It is evident from the plots that there was scatter in the wave motion, particularly 
at large displacements. However, there are some effects of rotation apparent. The 
inception point of a wave was defined as being that meridional displacement at which the 
wave crest first appeared. If an average inception point is computed for each test con- 
dition, is normalized with respect to the model radius, and is plotted against the spin 

parameter N /N , figure 12 results. Here it can be seen that for higher rotation rates, 
s p 

the position where the crest of the waves is first observed occurs at larger distances 
from the stagnation point. The movement of this point to larger meridional X values is 
concerned with the role the deceleration force has in forming the waves. McConnell 
(ref. 8) states that "under heating and the action of the deceleration force, a small 
amount of melt accumulates on the forward portion of the body. " Therefore, if the 
deceleration force is negated by centrifugal force, the inception point of the waves should 
be affected. The fact that the inception point of the waves changed with rotation rate 
necessitated the specification of the front surface to be based on the stationary model, 
as was indicated in a previous discussion. 

Further inspection of the waves reveals that a period of time and distance is required 
for them to accelerate to the relatively large values of velocity they ultimately attain. 




iz 

Time, t, sec 

(b) Spin parameter, 0. 1. 
Figure 11. - Wave crest position as function of time for various spin parameters. 
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Figure 12. - Effect of rotation on wave inception. 
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Since the position of wave inception is delayed at larger rotations, the surface area over 
which large wave velocities are present is decreased. Therefore, because the "melt 
movement and the wave motion are fundamentally related" (ref. 8), the effectiveness of 
the air shearing force in driving the liquid layer should be reduced. 


Temperature Rise of Body Fixed Point 

The temperature data, as measured by the thermocouples in the models, are plotted 
dimensionally in figure 13 and nondimensionally in figure 14. The basic parabolic shape 
of these curves is in agreement with the experimental results of reference 8. 

Investigation of figure 14 shows that the nondimensionalization did not correlate that 
data quite as well as it had the stagnation point recession results (fig. 9). It appears 
that the data for the runs in which the spin parameter equaled 2. 02 and 4. 2 are displaced 
a constant amount on the time axis from the remainder of the data. This would seem to 
indicate that the initial conditions for these two tests were different than for the others. 

A possible explanation is that small temperature gradients were present within these 
models, thus delaying the temperature increase in the vicinity of the thermocouple. 

These gradients could have been the result of the models manufacture because the length 
of tim e the models were inserted in the liquid nitrogen was not closely regulated. How- 
ever, even if the data from these tests were displaced the proper amount to the left in 
the plot, there would still be some scatter. McConnell (ref. 8) had a comparable amount 
of scatter in his dimensionless temperature data so that the results of this work are not 
unusual in that respect. 
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Temperature, 




Time, t, sec 


(b) Spin parameters, 0.88, 2.02, and 4. 2. 

Figure 13. - Temperature rise of body fixed point at various spin rates. 
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The conclusion of primary importance to be drawn from the data concerns that time 
during which steady-state ablation occurred. This period was characterized by the 
stagnation point receding at a constant rate. From figure 9, this is seen to occur 
approximately at dimensionless times greater than 3. 5. From figure 14, it is apparent 
that during this time the slopes of the curves faired through the data are all about the 
same. Therefore, it must be concluded that the temperature rise of the thermocouple, 
and hence, the rate at which heat was transferred to the vicinity of the thermocouple, was 
not affected by rotation. 


CONCLUDING REMARKS 

The effects of rotation on the ablation of an axisymmetric body composed of a highly 
viscous material and subject to deceleration in an airstream were studied analytically 
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and experimentally. The tests were conducted so that a parameter representing the im- 
portance of rotation effects and which was obtained from a nondimensionalization of the 
governing equations was varied from 0 to 4. 2. The following results were noted: 

1. The recession of the stagnation points of the models was not affected by rotation. 

2. Ablation in the vicinity of the stagnation point appeared to be unaffected by 
rotation. 

3. The heat transfer to a body fixed point on the major axis of the models was not 
affected by rotation. 

4. Liquid was propelled from the models in the vicinity of the flange when the spin 
parameter (N_/N_) equaled 0.88 or greater. 

5. From observation of data film, it was apparent that the amplitude of the waves 
decreased as rotation rate increased. 

6. The inception point of the waves tended to occur at larger distances from the 
stagnation point so that the surface area over which relatively large wave velocities 
occurred was smaller as the rotation rate increased. 

7. Rotation caused the waves to become more axisymmetric in nature. 

From these conclusions, it is apparent that the net effect of rotation on the ablation 
process for small meridional displacements was insignificant. This is understandable 
because of the smallness of the centrifugal force for these displacements. The temper- 
ature rise of the thermocouples also indicated that, over the lower portion of the forward 
surface of the models, rotation had no effect. At larger distances from the stagnation 
point, significant rotation effects on the characteristics of the waves and on the structure 
of the flange were apparent. 

Application of these results to the tektite problem indicates that the insensitivity of 
the ablation process to rotation over the front surface supports Chapman's arguments 
concerning entry velocities and angles. For rotation to have an effect away from the 
stagnation region on an actual tektite, the spin parameter (N g /Np) is approximately 1 or 
the rotation rate equals 2160 rpm, a seemingly prohibitive rate. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 29, 1967, 

124-09-03-01-22. 
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APPENDIX A 


SYMBOLS 


A acceleration 

Cp specific heat 

E modified Eckert number, 

D ?ef/< c p\< T i - 

e thermal layer thickness 

centrifugal force 

F component of centrifugal force 

cox 

in meridional direction 

F^y component of centrifugal force 

normal to liquid-gas interface 

g dimensionless acceleration 

parameter, p^AL/PgU^ 

H interface wave function 

k thermal conductivity 

L length characteristic of model 

size 

N index in viscosity -temperature 

relation 

N dimensionless pressure param- 

eter, P g U?.sfL/7T u U ref 

N modified Taylor number, 

^ 2 L% 2 /^H ref 

P pressure 

Pr Prandtl number, pC / K 

r 

p dimensionless pressure 

Q heat -transfer rate 

R radius of axisymmetric body 

Re Reynolds number, pUL/ p 


Re 0 rotation Reynolds number, 

P^ L /n u 

i dimensionless radius of axisym- 

metric body 

T temperature 

t dimensionless time 

t time 

U velocity in X-direction 

u dimensionless velocity in X- 

direction 

V velocity in Y-direction 

v dimensionless velocity inY- 

direction 

W dimensionless velocity in Z- 

direction 

W velocity in Z -direction 

X direction along liquid-vapor 

interface meridionally 

x dimensionless direction along 

liquid-vapor interface meridion- 
ally 

Y direction perpendicular to liquid- 

vapor interface 

y dimensionless direction perpen- 

dicular to liquid-vapor interface 

Z direction along liquid-vapor 

interface azimuthally 

z dimensionless direction along 

liquid-vapor interface 
azimuthally 
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rotation rate 


P 

7 

6 

V 

e 

M 

M 

P 

a 

T 


o 

convection parameter Re^ Pr^ 5^ 

ablation parameter 

scaling parameter for direction 
perpendicular to liquid-vapor 
interface 

dimensionless interface wave 
function 

dimensionless temperature 
dimensionless viscosity 
viscosity 
density 

time characteristic value 
shear stress 


w 

Subscripts: 

g gas 

I wave inception 

i interface 

l liquid 

m melting point 

o stagnation point 

ref reference 

T thermal layer 

oo center of model 

-oo gas free stream 
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APPENDIX B 


SUMMARY OF CHEN'S ANALYSIS 


Continuity Relation 


The continuity equation in dimensionless form is 


(U £ ) — (ru) + ( AlUref ) — (rv) = 0 
rel ax \ 5 1 I ay 


Since continuity requires that these two terms be of the same order of magnitude 


A 1 " 5 l 


Wave Motion Relation 


The wave motion at the liquid-gas interface is described by 


9H + U .M-V. =0 
at 1 ax 1 


where 


Y = H(x, t ) 


In nondimensional form, equation (B3) becomes 


+ 12- v, =0 


aU ref / 9t 1 9x 1 


where 


ii i ii i ii n 



I 


Since a requirement for wave motion is that drj/ 3t must be of the same order as the 
remaining terms, 

L - 1 (B7) 

* U ref 


Heat-Transfer Boundary Condition 


From a consideration of the heat-transfer rate across the liquid-gas interface, 


Qj - ” 

1 g 9Y 



3T 

9Y 


l; 


Nondimensionalizing results in 


K (T - T.) 

g v -oo i' 

de 


\ ( T i - t j~ 

de 

L e 

L g 

3y 'p 

g i 

l e i _ 

3y rp 


(B8) 


(B9) 


where 


T - T\ 

S T - T. 
-°o 1 


(BIO) 


and 


y T 


g L e a 


(Bll) 


since the two terms in equation (B9) must be of the same order of magnitude: 


K e ( T - 0 o- T i ) K J (T i- T =o> 


(B12) 
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Shear Stress Boundary Condition 


Equating the gas and liquid shear stress across the interface results in 


T i( x > *) = ^ ffi — 

1 S 1 9Y 


= M 


- au 


l i 


gi 


3Y 


(B13) 


l i 


The corresponding nondimensional equation is 


^■gi U -oo \ 9u 

L6 g Jay 


gi 


^ M Zi U ref \ 3u 
L6 7 Jay 


(B14) 


Z i 


The right and left sides of this equation are of the same order of magnitude so that 


ju .U u, .U f 
gi -°q ^ Z 1 ref 


(B15) 


Liquid-Boundary-Layer Considerations 

The energy equation (B14) provides a relation between the thermal and velocity layers 
in the body. In the steady state or long term, the amount of heat conducted into the body 
is of the same order of magnitude as that convected away in the moving layer. Therefore, 
it may be stated that 


(B16) 


Pr z Re z e z 


Chen provides a second relation between the thermal and velocity layers by assuming 
a linear temperature profile in the body or 


6, T. - T 

Z l m 


T i - T «, 

1 °° 


(B17) 
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Gas-Boundary-Layer Considerations 


The equations describing the heat-transfer and shear-stress relations between the 
liquid and gas at their interface inject two more unknowns into the problem; 5g and e^ 
These are determined from an order -of -magnitude analysis of the gas-boundary-layer 
equations. 

From the gas momentum equation, in order for the viscous terms to be of 1 order 
of magnitude, 


6 


g 



(B18) 


and from the gas energy equation 


e 


g 


1 

<Re g ) 1/2 Pr“ 


(B19) 


where 


m = 1 for Pr^ > 1 

m = — for Pr < 1 
2 g 


Solution for Characteristic Quantities 


The unknowns which remain in the problem are 6^, 
obtained by solving equations (B 7), (B12), (B15), (B16), 
follows: 


U re f> e i i and Tj. These are 
and (B17). The results are as 


Reference velocity: 


U 


ref 


“ U 


(C ) p r ( 1-m ) 

P; g 

(c p) 

g 



(B20) 
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Reference time: 


a = 


LPr Z < C p)„ 


U (Pr ) 
-°° v g 


(1-m) 


(CJ 


U 


T. - T 
1 

i T - T. 
m 1 


Scaling parameter for liquid velocity layer thickness: 


6 Z “ 


K, 




1/2 


T. - T 
j. m 

T - T. 

-oo x 


Scaling parameter for liquid thermal layer thickness: 


e l ~ 


K ; 

L K g (Pr g )m (R e g ) 1/2 _ 



Interface temperature: 


3 

y 


(ji rji3 

-oo " 1 i 


(T. - T )(T. - T ) 2 

v i oo /v i 


(B21) 


(B22) 


(B23) 


(B24) 


Ablation parameter: 


r 


r 

K f< c P >//"gi 

-\ 

1 

< 

K|(C p ) g P g Ma 

p r (3m-l) ” 

g 

J 


1/3 


(B25) 
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